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Abstract 

We study cosmological perturbations of self-accelerating universe solutions in the recently pro- 
posed nonlinear theory of massive gravity, with general matter content. While the broken diffeo- 
morphism invariance implies that there generically are 2 tensor, 2 vector and 2 scalar degrees of 
freedom in the gravity sector, we find that the scalar and vector degrees have vanishing kinetic 
terms and nonzero mass terms. Depending on their nonlinear behavior, this indicates either non- 
QJ dynamical nature of these degrees or strong couplings. Assuming the former, we integrate out 

• the 2 vector and 2 scalar degrees of freedom. We then find that in the scalar and vector sectors, 

gauge-invariant variables constructed from metric and matter perturbations have exactly the same 
quadratic action as in general relativity. The difference from general relativity arises only in the 
tensor sector, where the graviton mass modifies the dispersion relation of gravitational waves, with 
a time-dependent effective mass. This may lead to modification of stochastic gravitational wave 
spectrum. 
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I. INTRODUCTION 



The stability of general relativity (GR) predictions against small graviton mass has been 
a persistent challenge of classical field theory. The simplest ghost-free extension of GR with 
a linear mass term jj suffers from the van Dam-Veltman-Zakharov discontinuity, giving rise 
to different predictions for the classical tests in the vanishing mass limit {2, 3|. Although this 
problem can be alleviated by nonlinear terms the cost is the emergence of the Boulware- 
Deser (BD) ghost j^, which is generically unavoidable due to six degrees of freedom in the 
metric instead of the five of the massive spin 2 field. 

Adopting an effective field theory approach in the decouphng Umit^ the source of all the 
above issues can be traced back to the helicity mode of the graviton [6]. In this perspective, 
an analogue of the cancellation of the BD ghost in the linear theory by a specific choice of 
the mass term can be performed, rendering the nonlinear theory ghost-free up to quartic 
order by tuning the coefficients j'il]. 

More recently, a two parameter theory of nonlinear massive gravity was developed js, 9|. 
In this construction, terms to each order are chosen to remove the additional degree of 
freedom (would-be BD ghost) at the decoupling limit and it has the potential to be free of 



the BD ghost at fully nonlinear level 10l-ll2| . The more general construction with dynamical 



auxiliary metrics was also claimed to be free from the BD ghost at fully nonlinear order 
[l3|,Q. 



The massive extensions of GR are known to allow self-accelerating solutions |15l-l29|. 
This provides an application opportunity for the recent formulations of massive gravity, as 
an alternative approach to account for the current accelerated expansion. On the other 
hand, to model the accelerating universe in the scope of this infrared modified gravity 
theory, a cosmological solution is necessary. In the construction of with Minkowski 

fiducial metric, fiat Friedmann- Robertson- Walker (FRW) universe cannot be realized 22| . 
although cosmological solutions with negative spatial curvature exist {23! . For more general 
constructions with a nondynamical 30|] and dynamical 27|, |28| auxiliary metrics, maximally 
symmetric FRW with either curvature may be allowed. 

Given the nontrivial cosmological solutions which self-accelerate, it is thus a necessity to 
understand their properties against perturbations. In the light of the construction which 
removes the additional degree (would-be BD ghost), the expectation is to have five degrees 
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of freedom associated with the five polarizations of the massive graviton. On the other 
hand, the absence of BD ghost does not guarantee that the theory is safe; one still needs 
to determine the conditions under which the helicity and 1 modes, which are pure gauge 



in the massless theory, are stable 



3l|. Furthermore, these additional degrees should not be 



in conflict with observations. For instance, because of the emergence of a new scalar degree 
in the gravity sector, the Newtonian potential may acquire modifications from couplings 
between the matter sector and the helicity graviton. If this is the case, the parameters of 
the theory can be restricted by e.g. the solar system tests. 

The primary goal of the present paper is to address these questions. We consider the 
cosmological solutions in massive gravity as background and in the presence of a generic 
matter content, we present a gauge invariant formulation of perturbations. Although as we 
argued above, we expect the gravity sector to contain 5 dynamical degrees of freedom, at the 
level of the quadratic action, we show that the helicity and 1 graviton modes have vanishing 
kinetic terms but finite masses.^ As a result, in addition to the matter perturbations, the 
only dynamical degrees of freedom are the two tensor polarizations of gravity waves. 

This is exactly the same number of degrees of freedom as in GR. For the action 
quadratic in perturbations around the cosmological backgrounds, the mass terms of the 
potentially ghost-free construction turn out to be completely decoupled from the standard 
part (Einstein-Hilbert and matter terms) for scalar and vector modes. Thus, these modes 
evolve identically to their counterparts in GR, except for the additional cosmological con- 
stant contributed by the graviton mass. The only non trivial effect on the dynamics occurs 
in the tensor modes, which acquires a time dependent mass term determined by the fiducial 
metric of the theory. 

The paper is organized as follows. In Section II, we review the setup with a fiducial metric 
that is Minkowski and summarize the only cosmological solution 23|] allowed in this case. 
In Section III, we give a detailed study of the complete quadratic action for perturbations 
in the presence of a general fiducial metric of the FRW form, with an ansatz for the physical 
background metric also of FRW type with arbitrary spatial curvature. As an example, we 
specify to perturbations around the open universe solution [23] driven by a single scalar 
field matter in Section IV. We conclude with Section V, where we summarize our results. 



^ A similar situation was noted in the decoupling limit [l^, 0| • 
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The paper is supplemented by a number of Appendices, where the details of calculations are 
presented. 



II. OPEN FRW SOLUTION WITH MINKOWSKI FIDUCIAL METRIC 

In this section, we review the open FRW universe solution ||23il in nonlinear massive 
gravity [9| coupled to general matter content. 

The covariant action for the gravity sector is constructed out of the four dimensional 
metric gf^^ and the four scalar fields 0" (a = 0, 1, 2, 3) called Stiickelberg fields. The action 
respects the Poincare symmetry in the field space, i.e. invariance under the constant shift 
of each of 0° and the Lorentz transformation mixing them: 

(f)" ^(f)^ + 0" ^ A^0^ (1) 

The following line element in the field space is invariant under these transformations. 

Vabd(l)''d(j)'' = -{d(l)^f + 5ijd(j)'d(j)^. (2) 

Indeed, this is the unique geometrical quantity in the field space of 0". Thus the action can 
depend on (p"' only through the spacetime tensor 

= Vabd^rdA'. (3) 

In this language, general covariance is spontaneously broken by the vacuum expectation 
value (vev) of f^^. By assumption, matter fields propagate on the physical metric g^^, but 
are not coupled to f,^iu directly. The tensor f^,^, constructed from the invariant line element 
in the field space, is often called a fiducial metric. On the other hand, the spacetime metric 
g^j^iy, on which matter fields propagate, is often called a physical metric. 

The gravity action is the sum of the Einstein-Hilbert action (with the cosmological con- 
stant A) Ieh,a for the physical metric g^^, and the graviton mass term Imass specified below. 
Adding the matter action Imatter, the total action is 

I = Ieh,a [dfiv^cfi], (4) 

where 

WW = ^ j d'x^g{R-2K), (5) 

Imass[g|,uJ^lu\ = MpiTul j d'^ X^/^ {C2 + ot^C^ + aiCi) , (6) 
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and {(T/} (/ = 1, 2, ■ ■ ■) represent matter fields. Demanding the absence of gliost at least in 



the decoupling limit jsl, each contribution in the mass term Imass is constructed as 

A = ^([/Cf-[/C^]), 

£3 = i([/Cf-3[/C] [}C']+2[]C']), 

A = ^ {[JCf -6[ICf [IC'] + 3 [JC'] ' + 8 [/C] [JC'] - 6 [JC'] ) , (7) 
where the square brackets denote trace operation and 

K = si:-{VF^y^. (8) 

The square-root in this expression is the positive definite matrix defined through 

(VFv)'^ (v^)', = g^'fpu). (9) 

As already stated above, a vev of the tensor /^,^ breaks general covariance spontaneously. 
Thus, in order to find FRW cosmological solutions in this theory, we should adopt an ansatz 



in which not only (7^1, but also f^^ respects the symmetry of the FRW universes [23|]. Since 
the tensor /^i, is the puUback of the Minkowski metric in the field space to the physical 
spacetime, construction of such an ansatz is equivalent to finding a flat, closed, or open FRW 
coordinate system for the Minkowski line element. It is well known that the Minkowski line 
element does not admit a closed FRW chart but allows an open FRW chart. For this reason, 



in order to find open FRW solutions j23|], we first perform the field redefinition from 0" to 
new fields if"' so that f^^, written in terms of if"" manifestly has the symmetry of open FRW 
universes as 

= -n^iip')d^ip'd,ip' + a\ip')n,,iip'')d^ip'd,^\ (10) 

where i,j = 1, 2, 3, and 

QiAip ) = On H , — (11) 

is the metric of the maximally symmetric space with the curvature constant K (< 0). 
Concretely, this is achieved by 

0° = /(y.°)Vl-i^5.,W, 0^ = y=K/(y,V, (12) 

and 



n(v.°) = |/(/)|, a(^°) = v^|/(^°)|, (13) 



where / is a function to be determined and / represents its derivative. We then adopt the 
"unitary gauge" 



/ = t, (f' = x' 



so that 



ff^^dx^'dx" = -{f{t)fdi' + \K\ {f{t)ynij{x'')dx'dx^ . 



(14) 



(15) 



This is nothing but the Minkowski hne element in the open chart. For the physical metric, 
we adopt the open FRW ansatz 



(16) 



ds' = -Nitydr + aity nij{x'')dx'dx^ . 

Hereafter, we assume that > and a > 0, without loss of generality. 
The background action now yields, up to boundary terms, 

I = Mli j dt d^xNa^Vn {Leh[N, a] + m^L^,,, [A^, a, /]) + Irr^atter [N, a, a/] , (17) 



consisting of the Einstein-Hilbert part 

Leh = 

and the contribution from the mass term 



3K 3d2 



a2 a2 iV2 



(18) 



1 



6 + 4 ^3 + ^4 [S + b as + 2 a4) — (as + 



-sgnif/f) 



\f\a 



Na 

3 (3 + 3 ^3 + 0:4) 



(1 + 2 ^3 + 04) — [as + 04) 



(19) 



Hereafter, an overdot represents derivative w.r.t. the time t. 

Varying the action (ITTl) with respect to / yields the following constraint 



H-sgnU/f)- 



a 



X 



3 + 3 03 + 04 



( 1 + 2 0:3 + 04) — (0:3 + 0:4) 



a a 
where the Hubble expansion rate of the physical metric is defined as 

a 



0, (20) 



H 



Na 



(21) 
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Out of the three solutions of the constraint (1201) . the trivial solution a = sgn{f / f)\/—K N 
corresponds to the Minkowski spacetime in open chart. The remaining two branches of 



solutions are given by 



23| 



ait) = XMt) , ^ l + 2a, + a,±^l + a, + al-a, ^^2) 

as + a4 

and describe FRW cosmologies with K < 0.^ In the present paper we will focus only on 
these nontrivial cosmological solutions. 

Using the above constraint and varying the action (fT7|) with respect to N and a, we 
obtain the remaining background equations 

2H 2K 1 , , 

+ -r = T7riP + P)^ (23) 



N a? Mil 

where p and P are the energy density and the pressure of matter fields calculated from 

I matter 1 and 
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A+- ' 



[I + as) (2 + as + 2 - 3 a^) ±2 (l + as + a^ - a^f'^ . (24) 



(as + a^f 

Thus, for the cosmological solutions ( l22l) . the contribution from the graviton mass term 
Imass at the background level mimics a cosmological constant with the value A±. 

For a4 = (3 + 2as + 3a|)/4 and ±(1 + as) > 0, the effective cosmological constant K± 
vanishes, and the background solution reduces to the open FRW universe solution of GR. 
On the other hand, both X± and K± diverge for a4 = —as and ±(1 + as) > 0. In Figure [1], 
we show the sign of K± in the (as, a4) space. Note that X± are restricted to be positive by 
definition, as explained in footnote O Except for the restriction due to the positivity of X±, 



these are in agreement with the analogous region plots presented in Ref. 2J] ^. 



^ Note that X± are positive by definition since a{t) > and we assumed a{t) > 0. If we instead assumed 
a{t) < then the corresponding sohitions would be a(t) = —X±a{t) with the same X± and we would 
conclude X± > again. The essential reason for the positivity of X± is that the square-root in ([8]) is the 
positive one. 

^ Substituting — ^ 3 as, 12 a^, and switching the positive and negative branch definitions, our 

expression (p4)) recovers Eq.(6.6) of Ref. j23| . However, note that f^^, in the solution of 24] does not 
respect the FRW symmetry. 
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FIG. 1: Sign of the effective cosmological constant A± in the positive (left panel) and negative 
(right panel) branches. In the red (green) region with +45° (—45°) lines, Aj- is positive (negative). 
The white region and the dotted squared region correspond to 1 + 03 + a| — 04 < and X± < 0, 
respectively, and are excluded since the cosmological solutions (p2]) do not exist there. Along the 
dotted black line (defining the boundary between the red and green regions), A± = and the 
background solution reduces to the GR one. The solid line corresponds to X± = and thus defines 
one of the boundaries between the allowed (red or green) and excluded (dotted squared) regions. 
Along the dashed line, both X± and A± diverge, and it defines another boundary between the 
allowed (red or green) and excluded (dotted squared) regions. 

III. PERTURBATIONS IN GENERAL SETUP 

In this section we consider the graviton mass term Imassldfiu, f^u] defined by (I6])-(|8]) and 
( fTO|l -f lTT|) . but with an arbitrary value of K and arbitrary functions n{ip^) and a(v?°)- We 
shall develop a formalism to analyze perturbations of this generalized system around flat 
[K = 0), closed {K > 0) and open {K < 0) FRW universes. Cosmological implications of 
this type of generalized massive gravity will be discussed in future publication. 

For the background we adopt the physical metric of the FRW form f[TB|) with f lT^ . but 
with general K, n{ip^) and Q;(v?°). Without loss of generality, we assume that > 0, n > 0, 
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a > and a > at least in the vicinity of the time of interest, where and a are the 
lapse function and the scale factor of the background FRW physical metric. (Otherwise, 
we consider |A^|, \n\, \a\ and |a| and rename them as N, n, a and a.) As reviewed in 
the previous section for open universes in the case of the Minkowski fiducial metric and as 
shown in Appendix |X] for general cases with arbitrary K, n{ip^) and Q;(v5°), the background 
equation of motion for the Stiickelberg fields (p"" has three branches of solutions. One of 
them does not allow nontrivial cosmologies and thus is not of our interest. The other two 
branches of solutions allow nontrivial cosmologies and are given by (!22|) even for general K, 
n{ip^) and Q;(v?'')- In this section we then consider perturbations of the physical metric and 
the Stiickelberg fields around the FRW solutions in these nontrivial branches. 

A. Exponential map and Lie derivative 

Since the fiducial metric /^;^ is defined as in ( fTOll without referring to the physical metric 
g^y, let us begin with perturbations of the Stiickelberg fields ip"". We define perturbations 
vr" of if"" through the so-called exponential map. Actually, since the action will be expanded 
only up to the quadratic order, we can truncate the exponential map at the second order. 
We thus define vr" by 

(^» = a;" + 7r" + i7r^a;,7r" + 0(e^), (25) 

or equivalently, 

r,- = (^« - x") _ !(/ - x')d,{^'' - x") + 0{e'). (26) 

Here, e is a small number counting the order of perturbative expansion: tt" = 0(e) and 
_ _ (9(^g) By substituting the expansion (!25|) to the definition of the fiducial metric 

Uu = fab{v')d,'p''d,ip\ (27) 

where 

/oo(^^) = -^2(^0), /o,(^^) = ^(^^) = 0, U,{^^)=a\^')n,,{^''), (28) 

we obtain 

f,u = Uu{x') + Cj,,{xn + \ (^.)' UAx") + 0(e3). (29) 

Here, C^^ represents the Lie derivative along vr^. Actually, this formula is not restricted to 
(EHD but holds for any fab{^^)- 
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B. Stiickelberg fields and gauge invariant variables 



We define perturbations 0, /3i and hij of tfie pliysical metric by 

goo = -iV2(t)[l + 20], 

goi = N{t)a{t)(5i 

gij = a^(t) [Qijix'') + hij] . 

We suppose tliat 0, A, hij = 0(e). 

Under tlie linear gauge transformation 

^ x'^ + C^, iC" = 0(e)) 

each variable transforms as 



TTj TTj + ^i, 

^ 0+^9i(iVe 



A ^ A - - + l^e., 

a i\ 
hij h,j + D,^j + Dj^i + 2NH^'^Qij, 

where H is the Hubble expansion rate as defined in (|2T1) . 

and is the spatial covariant derivative compatible with Qij. 
We then define gauge invariant variables 

r = 0-^5,(iVvr°), 
a A/ 

For later convenience, let us decompose /3f and h^j as 

= D,p- + SJ, 
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where and F[ are transverse, and 'jtj is transverse and traceless: 

D'S^ = D'FJ' = 0, DSj = 0, n'H, = 0, (36) 



and D' = Vl'Wj. 



C. Graviton mass term 

At the FRW background level, the graviton mass term acts as an effective cosmological 
constant A± shown in (p^ . The proof of this statement is presented in Appendix Rl for 
arbitrary K, n{ip^) and a Thus, calculations are expected to be simplified if we add 
Mpi J d'^x^y—gA± to Imass before performing perturbative expansion. For this reason, we 
define 

^massld ^J,Iy 1 ffJ-u] = I mass[9fJ.u^ ffiu] ~^ M PI J d'^Xy/—gA±, (37) 

and expand it instead of Imass itself. 

As shown explicitly in Appendix |Al upon using the background equation of motion for 
the Stiickelberg fields but without using the background equation of motion for the physical 
metric, the graviton mass term can be expanded up to the quadratic order as 

Imass = Imass + ^massi^ij] ~^ )) 

I^Lm = ^ / d'xNa'VnM'aw [{hn' - Khl] , (38) 
where the zero-th order part Imass is independent of perturbations. 



= ±{r - l)mlXl^Jl + ^3 + ai - 04, (39) 

na 1 H a a 

r = = , H = , Hf = — , 

Na X±Hf' Na ^ na 



X± is given by (1221) . and 

''IJ 5 ' ''IT 

With the decomposition of h^, in (!35|) . the quadratic mass term is expanded as 



/i" = Q'^h^., = Q^'^Q^^hli. (40) 



41s = Mj,iJ d'xNa'VnM, 



GW 



3{ry - ^E-^iA + 3K)E- + ^F;(A + 2K)F^ - ^Y^Hj 



(41) 
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where 

= n^^Fj, = n^'n^'^ki. (42) 

What is important here is that the quadratic part ImLs is gauge- invariant and depends 
only on hj-, or equivalently E'^ , F[ , 'jij). In particular, it does not contribute to the 
equations of motion for and 

We note that Mq^^ vanishes or diverges for some special values of the parameters (as, 
a4): 

a4 = -3 (1 + as), ±(as + 2) > ^ M^^ = , 

^4 = 1 + as + a^ =^ Ml^ = , (43) 
a4 -as, ±(l + as)>0 =^ |M^^y| oo, 

where the ± signs are for the ± branches, respectively. In the following we suppose that the 
parameters (as, a4) take generic values away from the special values shown in (H3ll . 



D. Matter perturbations and gauge- invariant variables 



Let us divide matter fields o"/ (/ = 1, 2, ■ ■ ■ ) into the background values af^ and pertur- 
bations as 

aj = aP + 5aj. (44) 

We suppose that {cr/} forms a set of mutually independent physical degrees of freedom. 
Otherwise, we consider a subset of the original {cr/} consisting of independent physical 
degrees of freedom and rename it as {cr/}. We can construct gauge-invariant variables Qi 
from 5ai and metric perturbations, without referring to the Stiickelberg fields. 
For illustrative purpose let us decompose hij and as 

Kj = 2^a, + (^DiD, - ^n^JA^ E + ^{D,Fj + DjFi) + 7,,-, 

= Ae + eo (45) 

where Si, Fi and are transverse, and A is the Laplacian associated with Vtif 

D'Si = D'Fi = D'iJ = 0, A = D'Di. (46) 
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Under the gauge transformation (13T!) . each component of the physical metric perturbation 
transforms as 



a i\ 

E E + 
Si ^ Si + 
^ + 2ef , 



(47) 



Noting that the vector defined by 



(48) 



transforms as 



Zf" + 



(49) 



we can construct the following gauge-invariant variables out of matter perturbations and 
physical metric perturbations: 



Qi = 6a J - Cz(yf\ 
$ ^ 0-l9,(iVZO), 

^> = i)- NHZ° - -AE, 
6 

Bi = Si —Fi , 

2N 



(50) 



and jij is gauge- invariant by itself. In the above, Cz is the Lie derivative along Z^^. 

Those gauge-invariant variables defined here and in Subsection IIII B\ i.e. {Qi, ^E', Bi, 
jij, (j)'^, (3'", S'f , ip'", E'^, F^}, are not independent. Indeed, it is easy to show that 



-— ( V" - ^ - -AeA + —E\ 
aHV 6 2N ' 



SI = B. + -Fr 



(51) 
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There are no more independent relations among gauge-invariant variables defined here and 
in Subsection IIIIBi ^ Therefore, we have the following set of independent gauge-invariant 
variables. 

{g,, $, ^, Bi, 7,,, r, E\ F^. (52) 

Based on their origins, we can divide this set of independent gauge-invariant variables into 
two categories as 

{Qj, $, v]/, B,, 7,,} and {r, E\ F^}- (53) 

The first category consists of those gauge-invariant variables that originate from the physical 
metric g^^ and the matter fields {o"/}. Thus, those in the first category already exist in GR 
coupled to the same matter content. On the other hand, those in the second category are 
physical degrees of freedom associated with the four Stiickelberg fields ip"". 

E. Structure of total quadratic action 

Let us now define 

i[9^lu^ CT/] = IEH:A\9^lu] + Iraatter[g^lu, CT/] , A = A + A±, (54) 

so that 

I = I[9fiiy, CT/] + Imass[gtiu, f'liu]- (55) 

Since Imass was already shown to be gauge-invariant up to the quadratic order, (l55|l implies 
that I is also gauge-invariant up to that order. Thus the quadratic part F"^^ of / can be 
written in terms of gauge-invariant variables constructed solely from perturbations of the 
physical metric perturbations (0, hy) and matter perturbations 6aj, i.e. {Qi, Bi, 

Therefore, the total quadratic action has the following structure. 

= /(2)[Q,,<|.,vI/,i?,,7d E\ i^^7d, (56) 

where the explicit form of ImLs is shown in fj4T]) . As already stated, gauge-invariant variables 
listed in f l52|) are independent from each other. 

^ See also the sentence just after (H^ . 
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Note that ip'^, and FJ^ do not have kinetic terms but have nonvanishing masses, 
provided that the parameters (03, 04) take generic values away from the special values 
shown in Thus, we can integrate them out: their equations of motion lead to 

= = 0, FJ" = 0, (57) 

and then 

/If ^ r 

J(2) = J(2) [Qj, $, vl>, 5„ 7,^.] - ^ y d'xNa'VnM'awl'^^.j, (58) 

where Mq^t is given by f l5I?]) . For scalar and vector modes, this quadratic action is exactly 
the same as that in GR with the matter content {cr/}. 



F. Gravitational waves with time-dependent mass 

The total quadratic action for the tensor sector is 

1 1 



(59) 



provided that there is no tensor- type contribution from the quadratic part of Imatter- In 
this way the dispersion relation of gravitational waves is modified. The squared mass of 
gravitational waves M^y^ is given by ( 139|) and is time-dependent. 

If MQ^r is negative then long wavelength gravity waves exhibit linear instability. For 
generic values of parameters (0:3, 04) away from the special values shown in (H3l) . we see 
from the formula fl5^ that the sign of M^^y is the same as the sign of the combination 
±(r — l)mg, where ± signs correspond to ± branches, respectively. 



IV. AN EXAMPLE: SCALAR MATTER FIELD AND MINKOWSKI FIDUCIAL 

In the previous section we have analyzed quadratic action for perturbations around non- 
trivial FRW backgrounds, with a general FRW fiducial metric and a general matter content. 
For scalar and vector modes, we have shown that the quadratic action is exactly the same 
as that in GR with the matter content. For tensor modes, on the other hand, we have seen 
that gravitational waves obtain a time-dependent mass. 

In this section, in order to illustrate these results a little more explicitly, we shall consider 
a simple example consisting of the massive gravity with Minkowski fiducial metric, coupled 
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to a canonical scalar matter field with potential V. The total action of this system is 

'R o,. . .A 1 



Mj>i ( — + mg(£2 + asA + a^Ci] 



(60) 



where Ci^2,z are the graviton mass terms defined in ([7]). In the case of the Minkowski fiducial, 



the only nontrivial FRW background is the open FRW solution found in [23|] and reviewed 
in Sec. [Tll Thus, in this section the curvature constant K is set to be negative and the form 
of the fiducial metric is specified by (|T0l)-( |T3l) . 

For the FRW background (fT^ - (fTB|l with K < and a = cr'^°^(t), the equations of motion 
read 



3H^ + 



3K 



A± + 



(0)\2 



2iV2 



2H 2K 

+ -IT 



(0)\2 



Mil m 



1^ 



0. 



(61) 



We now introduce perturbations to the metric g^^, the four Stiickelberg fields cp"" and the 
scalar matter field a. We will be developing the perturbation theory without specifying a 
gauge and in the end we will switch to gauge invariant perturbations, as we have already 
done in the previous section in a more general setup. The total quadratic action before 
switching to gauge invariant perturbations is presented in Appendix [Bl We then adopt the 
decomposition of the form 



TTj = DiTi + vrf , 



A = + Si 



(62) 



where ttJ , Si and Fi are transverse: 



D'-kJ = D'Si = D'Fi = , 



(63) 



A. Tensor sector 



We start by considering the tensor sector. We use the decomposition fl62p in the total 
action (IBl|) - flB3p . keeping only the transverse traceless mode 'jij. Since the matter sector 
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has no tensor degrees, the action is the same as the one given in (159|) with (139|) and for the 
Minkowski fiducial, we have 

Notice that in the accelerating universe, r and thus Mq^^ grow at late time. 
Switching to conformal time drj = N dt/a and defining the canonical fields 

lij = lij , (65) 



the tensor action takes the form 



(66) 



where a prime denotes differentiation with respect to conformal time. Next, we use harmonic 
expansion through 

7ij = / k'^dk'y^Yij{k,x) , (67) 



where Yij{k,x) is the tensor harmonic satisfying 

(^A + k^^ Yij = , D%j = 0, n'^Yij = 0, (P = %fP) (68) 

with k"^ > \K\ taking continuous values. Suppressing the momentum index, we obtain the 
equation of motion 



a" 



j" + \^k^ - — + 2K + a'M^wjj = 0. (69) 

As we showed in the previous section for a generic setup, the tensor mode acquires a mass 
contribution whose time dependence is determined by the fiducial metric (which is Minkowski 
in the present example). 



B. Vector sector 

We now move on to the vector sector. Keeping only the transverse vector modes in the 
decomposition (l62i) . the total action ( ]Bll) -( lB3|) reduces to 

^vector ~ g^' J d^xNa^y/^i elector (70) 
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with 



elector = ^ [d'FW.F, -2KF'F^ - — [d'FW.S, -2KrS, 

M^aw {D'FW.F, - 2 KF^F,) + ^ D'SW.S, - ^ S,S' 
-2 M^^ ( Avrf + 2 K ttJ) [F^ - vr^ , (71) 



where 



S' = Q'^Sj, F' = Q'^Fj, 4 = Q'^nJ. (72) 
We then switch to gauge invariant perturbations, as defined in (lM|) - (!35l) . and obtain 

elector = ^ [d^FID^FJ -2KFIF:) - ^ [d^FID.S] -2KFIS^) 

-]- M'aw {D'FiD.FJ - 2 KFIF-) + -| D^SiD.S^ - ^ S^SJ , (73) 



which is manifestly gauge invariant. Varying this action with respect to yields an alge- 
braic equation for S^, which can be solved by 

SI = ^^r (74) 
Using this solution back in the action, we get 

4eL = ^ / d'xNa'VnM'^^ F;(A + 2K)FJ^ . (75) 

This clearly shows that the kinetic term for vector perturbation vanishes at quadratic order. 
Provided that the parameters (03, ^4) take generic values away from the special values 
shown in ( l43l) . the equation of motion for F^ leads to F^ = and then ^Jctor = 0- 

C. Scalar sector 

Finally, we consider the scalar perturbations. Using the decomposition (1621) in the total 
action ( lBll) -( |B3l) . we obtain the action for the scalar sector as 

= ^ / d'xNa'VnCscalar , (76) 
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with 

^scalar 



_ 1 D.^^jji^^ _ V" - ^ (0 - 3 ^) 5(T - 2 l^' (0 + 3 i/j) 6a 



K 



12 H 

N 



12 K 



K 



6 
2K 

2 



K 



Mil m 



Dil3 D'I3 + 



aN 



a 



2/\ip-- A'E ~KAE-2 NHAcp + A6a ] (3 
3 Mpi 



2 All)- - - K + 6 [Ml 



K 



+ — 2 D(i\)D'i) + - A^AE + — DiAE D'AE -2KDi^ D'E 
\ 3 18 

+2 M^^ (^K ( Avr D'n - DiED'n) + ^AE An - 24^ An - 6 H N n'^^ 



+2 H Mf 



GW 



N 



Avr + 3 Hn^ n° 



(77) 



Hereafter in this subsection, a represents the background value a^^\t). Next, we switch to 
gauge invariant variables defined in ( |Ml) - (l35ll and carry out harmonic expansion. We then 
obtain the equations of motion for the nondynamical degrees as 



1 f K k'^ -SK ■ 



-3H' + 



6N 



0" 



2Mj,iN 



6a^ 



2 Mi^m 



,^ k'^-?>K 3H . k^(e-3K)^^ 



N 



6 a' 



2 Ml, 



2 MliN'^ 



(78) 



where 



= 6a- an\ 



(79) 



and substitute their solutions into the action. Defining the analogue of the Sasaki-Mukhanov 
variable 



a 



NH 



(80) 



the resulting quadratic action reads, up to boundary terms. 



y(2) 

scalar 



Skdta^ N 



6 Miw\r? - g Mi^e [e - 3 k) i^^^p + Cq 



(81) 
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where 



and 



2H^P-3K) ( 1 ,-2 



83) 



= y" _| i \ _| Lj^ L : 

« a^V KMl,m) KMl,H(2H^{k^-?,K) + K^) 

Provided that the parameters (as, take generic values away from the special values 
shown in the equation of motion for ip'^ and E'^ lead to = E'" = and then 

lilar = ^ld'kdta'NCc. (84) 

As shown in Appendix [Cl the action (18^ agrees with the standard results in GR coupled 
to the same scalar matter field a. 

To summarize, the scalar sector consists of a dynamical degree which evolves exactly like 
the standard Sasaki-Mukhanov variable in GR, and two more degrees which have infinite 
mass. 



V. SUMMARY AND DISCUSSIONS 

In the context of the potentially ghost-free, nonlinear massive gravity and for a gen- 
eral fiducial metric of FRW form, we analyzed linear perturbations around self-accelerating 
cosmological solutions of arbitrary spatial curvature, populated by generic matter content 
that is minimally coupled to gravity. By constructing a gauge invariant formulation of per- 
turbations, we found that massive graviton modes in the scalar and vector sectors have 
vanishing kinetic terms but nonzero mass terms. By integrating them out, we showed that 
the part of the action quadratic in scalar and vector perturbations is exactly the same as 
in GR with the same matter content. In other words, the dynamical degrees in the gravity 
sector comprise only the two gravity wave polarizations. We also found that these acquire 
a mass whose time dependence is set by the fiducial metric. 

In Fierz-Pauli theory in de Sitter background, it has been known that the scalar mode 
among five degrees of freedom of massive spin-2 graviton becomes ghost for 2H^ > m\p, 
where H is the Hubble expansion rate and mpp is the graviton mass 3l|. This conclusion 
does not hold in the nontrivial cosmological branches of the nonlinear massive gravity. In- 
deed, as stated above, the scalar and vector modes have vanishing kinetic terms and nonzero 
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mass terms for any FRW background. This sharp contrast to the hnear (Fierz-Pauh) mas- 
sive gravity stems from a pecuhar structure of the graviton mass term expanded up to the 
quadratic order in perturbations: it depends only on the (ij)-components of metric per- 
turbations and thus are independent of (00) and (Oz)-components. This Lorentz-violating 
structure is possible because the vev of f^i, in the cosmological branches spontaneously 
breaks diffeomorphism invariance in a nontrivial way. 

It is still fair to say that the nature of the cancellation of the kinetic terms in the quadratic 
action of scalar and vector sectors is not well understood. This may be an indication that 
these sectors exhibit an infinitely strong coupling. If this is the case then we cannot properly 
describe the scalar and vector sectors without knowledge of a UV completion. On the 
other hand, since the modes have non vanishing masses, it may be possible that these are 
just infinitely heavy modes without low-energy dynamics, as we have assumed in the main 
part of the present paper. In this case we can safely integrate out those extra modes and 
trust the resulting low energy effective theory. In order to establish the fate of those extra 
degrees of freedom, i.e to judge whether they are strongly-coupled or non-dynamical, the 
linear perturbation theory is not sufficient and nonlinear methods are needed. This study is 
beyond the scope of the present paper and is left for a future work. 

On the other hand, the modification in the tensor sector may leave a signature in the 
stochastic gravitational wave spectrum. The additional term in the mass of the tensor 
modes is time dependent, while its sign is determined by both the fiducial metric and the 
cosmological evolution. A positive but large contribution may give rise to a suppression 
of the gravity waves and a null signal in the large scale tensor-to-scalar ratio. However, 
this deviation from scale invariance may allow the signal at small scales to be potentially 
observable in the space-based gravity wave observatories such as DECIGO [32], BBO 33] 



and LISA 



Mi- 



lt is important to note that the present analysis is purely classical and special care 
is needed when discussing the evolution of cosmological perturbations which start off in 
quantum mechanical vacuum. In order to address quantum issues such as the radiative 
stability of the structure of the effective theory describing cosmological perturbations, one 
of the first important steps is to identify the strong coupling scale below which the nonlinear 
massive gravity theory can be trusted. In the branch described by the trivial solution to 
the Stiickelberg equation of motion (I2U]) . the so called decoupling limit has been useful for 
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this purpose. However, since this trivial branch is not compatible with FRW cosmologies, in 
the present paper we have considered other two branches described by non-trivial solutions 
In these cosmological branches, the usual decoupling limit is not applicable, at least 
apparently. We thus need to develop a new technique to identify the strong coupling scale, 
or directly analyze nonlinear dynamics of the whole system. This is certainly one of the 
most important issues in the future research. 
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Appendix A: Calculation of graviton mass term 

In this appendix we consider the graviton mass term introduced in Sec. Hi] and generalized 
at the beginning of Sec. Illli We thus consider arbitrary K, n{ip^) and a{(p^)- Without loss of 
generality, we assume that A^>0, r;,>0, a>0 and a > at least in the vicinity of the time 
of interest, where N and a are the lapse function and the scale factor for the background 
FRW physical metric. 

We expand the generalized graviton mass term up to the quadratic order in perturbations. 
In doing so, we shall not use the background equations of motion for the physical metric 
since they depend not only on the graviton mass term but also on the Einstein-Hilbert action 
and the matter action. On the other hand, we shall use the background equation of motion 
for the Stiickelberg fields in the middle of calculation. 

1. Background equation of motion for Stiickelberg fields 

In this subsection we shall derive the background equation of motion for the Stiickelberg 
fields 9?" by expanding the graviton mass term Imass up to the linear order in tt" (= 5^9") 
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without variation of the physical metric (yf^j^. This will be a good warm-up for the forthcoming 
subsections. 

Using the formula ( 12^ . /^,^ is expanded up to the linear order as 



/( 



00 



—n 



1 + -5t(n7r°) + 0(e' 



n 



foi = fio = otn 



a n 



f,j = a" [(1 + 2nHfTT'^)nij + Att^ + D,tt, + 0{e'')] 



where an overdot represents differentiation with respect to the time t and 

a 



na 



With the unperturbed physical metric 

5'oo = 



^^ 9oi = 9io = 0, gij = a^VLi 



this leads to the following expansion for /'^ (= g^^ ft 



pv) 



J 



n 



1 + -a(n7r°) + O(e^) 



an 



n 



a 



n 



J « ]\f2 



a 



n 



f. = — [(1 + 2nHf'Kyi + D'ttj + Djtt' + 0{t 



(Al) 



(A2) 



(A3) 



(A4) 



Then, using the formula in Appendix [D] for matrix square-root, /C^^, defined by (IE]) 
expanded up to the linear order as 



IS 



(A5) 



where 



and 



(0)0 



(0)0 



0, /C 



{0)i 



0, }C 



(0)i 



(1)0 



(1)0 



na 



n 



n r. a . 

Avr° + -VTi 

a n 

U TT H TT 



(l)i 



a(l + r) L a n 



(A6) 



(AT) 
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Here, we have defined 

It is now straightforward to expand the graviton mass term ([6]) up to the first order. The 
result is 

d^xNa'Vn—iaH - aHf)J^n'' + 0{e^), (A9) 
where the zero-th order part JmLs does not depend on vr'^ and 

J0 = 3-2X + 03(1 -X)(3-X) + 04(1 ^ = (AlO) 
Therefore, the background equation of motion for the Stiickelberg fields is 

{aH - aHf)J^ = 0, (All) 

where Hf is defined in f lA2p and 

H = —. (A12) 

Setting aH = aHf would not allow nontrivial cosmologies since in this case the background 
evolution of the physical metric would be determined not by the matter content but by 
the fiducial metric. Thus, we restrict our attention to solutions of = 0. This leads to 
X = X±, where X± are given by (l22l) . 



2. Unitary gauge 

In this subsection we shall expand the mass term up to the quadratic order in the unitary 
gauge, i.e. under the gauge condition 

7r° = TT^ = 0. (A13) 

After obtaining the expression in the unitary gauge, it is relatively easy to infer the corre- 
sponding expression in general gauge and to confirm it. The expression in general gauge 
will be presented in the next subsection. 

In the unitary gauge, the fiducial metric f^u is the same as in the background. Hence, 
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/'^ (= g^^fpv) is expanded up to the second order as 



fo 



n 

Na 



[l_20+(402_/3*/3.) + O(e3)] 
[/3._20/3,-/3^-/i,, + O(e=')], 



n 



Na 



a 



where 



(3' = n'^/3j 

Then, taking the matrix square-root as prescribed in Appendix [Dl K.^^^ defined by 
expanded up to the quadratic order as 



where 



~ 



(1)0 



n 



/C 



(1)0 



nr(3'^ 



r«'' 

N{l + ry ° a(l + r) 



^ 2a^' 



and 



r(2)o 

^(2)0^ 



iV 



3 . r 2 + r . ^ 
2^ 2 1 + r '^'^ 



a 



N{l+r) 
nr 
a(l + r) 



1 + r 2(1 + r 



(2)i 



a 



2 + r 
l + r^ 2(1 +r 

kj 



Pkh' 



(AM) 
(A15) 

i-dni) is 

(A16) 

(A17) 
(A18) 



(A19) 



_2(l + r)2'^ 

It is then straightforward to expand [/C"] (n = 1,2,3,4) up to the quadratic order. The 
result is 

[/C"] = [/C"](°) + [fC^f^ + [/C"](') + 0{e'), (A20) 

where 



l(0) 



[/C"]^^^ = nrX(l -rX)"-V+ ^X(l -X)"-^/i, 



3(1 -X)" + (1 -rX)^ 



(A21) 
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and 



[ICY^> = 


3 , 

— rJ 
2 


[ic'Y'^ = 


(4rX 


[K'f^ = 






3X 







^ij ~l~ o P Pi) 

2ri 
3 
4 



(r2-2r3X + r4X^)/3*A, 



X 
— ( 

ri 



41 (2) 



6(2rX - 1)(1 - rX)VX02 + _(2X - 1)(1 - XfXh'^h 



+ 



2X 



(r2 - 3r^X + 3r,X'' - r.X^)^'^^. 



ri 



Here, X is defined in ( lAlOl) and 



(A22) 



(A23) 



j=0 



We are now ready to expand the graviton mass term ([6]) up to the quadratic order. 
However, before doing so, let us expand it up to the hnear order in order to see the effective 
energy density due to the graviton mass term. The result is 

l.\ 1 



-{<l>+-h P, + -M'pim'il-r)XhJ, 



+ 0(e"), (A24) 



where the zero-th order part Imlss does not depend on perturbations, and 

Pg = -Ml,ml{l - X) [3(2 - X) + a3(l - X)(4 - X) + a,{l - X)^] 



(A25) 



is the effective energy density due to the graviton mass. 

Having obtained the expression for pg and noticing that the factor (0 + /i/2) in (1A24P is 
the linear order part of 



-9 



we expect that expanding 



+ 0(e=^), (A26) 



(A27) 



instead of Imass itself, should simplify the resulting expression. We thus expand Imass up to 
the quadratic order. The result is 



imass = Mliml I d'^xNa'^Vn 



mass ' mass ' mass ' ^ \ / 



(A28) 
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where 



^^^L = -^^(1 - ^) [3 + 3«3(1 -X) + a,il - Xf] 



L^±s = \{l-r)XhJ^, 



f (2) 
mass 



(t)h + — + -(1 - r){h' - 2h'^h 
1 + r 4^ 



1 



4. 

2 n /r2 / 7 2 



and 



+ -m;^M^;^(/i^-/i*^/i,,), (A29) 



m-^Mlw = XJ^ + (1 - r)X2 [1 + a^{2 - X) + a^{l - X)] . (A30) 



As shown in the previous subsection, the background equation of motion for the 
Stiickelberg fields for nontrivial cosmo logical branches is J,^ = and gives X = X^., where 
J0 and X± are defined in flAlOp and fl22l) . respectively. For X = X±, it is easy to see that 



Pg = Mp;A±, where A± are defined in (!24l) . Thus, upon using = 0, the definition of Imass 
in flA27p reduces to fl371) . Also, for X = X±, defined above reduces to that defined in 

In summary, in the unitary gauge, upon using the background equation of motion J,/, = 
for the Stiickelberg fields but without using the background equation of motion for the 
physical metric, Imass defined in (!37|) is expanded up to the quadratic order as 

/M'2 r 

Imass[g,u, f,u] = 41. + d'xNa' ^QM^awih' - h'^K,) + 0{e'), (A31) 

where the zero-th order part Imass = Mpim'^ J d'^xNa^y/TlLmass does not depend on the 
perturbations. 



3. General gauge 

In general gauge, we expect that the expansion of Imass defined in (1A27I) should be similar 
to that in the unitary gauge, provided that each metric perturbation variable is replaced 
by the corresponding gauge-invariant variable constructed from the metric perturbation and 
the Stiickelberg field perturbation. Such gauge- invariant variables are defined in (134|) . We 
thus expect that 

Imass[g,u, Uu] = 41. + dx' N a' V^M^^ [{h^f - h^^hl] + A^,,„ (A32) 
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with relatively simple expression for Amass up to the quadratic order. Of course, for J^j, = 
tt" = 0, Amass should vauish up to the quadratic order, as shown in the previous subsection. 
Indeed, by direct computation we can confirm that this is true not only in unitary gauge 
tt" = 0, but also for tt" 7^ as far as = is imposed: 



A = Ofe'^ 

'—^mass ^ \^ 



for Ja, = 0. 



(A33) 



If we do not impose = nor vr" = then A^ass up to the quadratic order is 

1 — r 



1 

+ 2 



2N 



'1 - r)(h + SNHhy - —(1 - r)/3'Di7i° 

a 



N^{l + r) 



+3(1 - r)7r" ( -^dtiNa^Hn^) + NHDiTr 



X 



+ - 



-(20 + h) + —-dt{Na\'^) + D.iTx' 



n'^X^Js + Oi^). 



(A34) 



Appendix B: Quadratic action for the scalar field example 



This appendix presents the total quadratic action for the scalar field example of Sec. IIVI 
Combining the graviton mass term ( ]A32[) -( !A33I) with the Einstein-Hilbert and scalar field 
parts, the complete quadratic action can be calculated up to boundary terms as 



= Mil I d^xNaWniC + 



-9 



(2) 



3H' + 



3K 



A 



± 



a 



Mil \2N^ 



V 



+ {h^-2hijh'^) ( _-_ + 



M|;A^2 



(Bl) 



where 



-9 



(2) 



')2 1 1 11 

- + -(3i(3' - -hiA'^ + -h'^ + -(l)h. 



(B2) 
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and 
C 



^ ( kA'^ - + E(i)h - - ( 2 H (f) - — h] Dil3' — Dil3. h'i -3H 



2 J,2 



4a2 

+ [2 D'hik DjW^ - Dkhij D^h'^ + 2h DiDj0 ~ h Ah] - ^ [hijh'^ + Ah 



8a2 

2 

GW 



4a2 

- (h^ - h,jh'^) - - (Din' + 2NH7r^)h + - h^j D'tt^ + 3 (11^ 
8 2 2 



[{DiTv'y - Di-KjD'-K^ + 2K TTiTT*] + 2HN n^Diir' 



Ml, L 



2a^ aN'^ 2 



V 



(B3) 



Appendix C: Comparison with GR 



In this appendix we compare the scalar action ( 184|) with the GR results in the literature. 
For this purpose it is useful to define the gauge invariant Bardeen potential, which in our 
language corresponds to 



^ N dt\ ' 2N 



(CI) 



or, using Eqs.( l6T]) . (1781) and ( IHOl) . as well as the equation of motion for Q obtained from 
varying fl82p . we can write 



a 



2 Mil m (A;2 - 3K) + Kj^ N 



Q 

N 



3H 



2 Mil H m 



+ I Q 



a 



(C2) 



Using the above definition, along with the equation of motion for Q, then switching to 
conformal time adrj = N dt, we see that the equation of motion for $ in the standard 
scenario, given in Eq.(8.140) of Ref. 35|, is satisfied. To compare the normalization of the 
action with the literature, it is also useful to define the following quantities 



2 a Ml 



PI 



NH^ KHN'^ 



(C3) 



H a 

With these definitions, it is straightforward to verify that the equations of motion given 
in Eqs. (21-22) of Ref. jsGj are satisfied, with substitution dt — )■ Ndt. Furthermore, using 
the above expressions, the action given in Eq. (24) of 3^ reproduces our action (l84l) up to 
boundary terms. 
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Appendix D: Perturbative expansion of matrix square root 

Let us expand an x matrix A as 

A = + eA(i) + eM(2) + Oie'), (Dl) 
where the zero-th order part is assumed to be of the form 

(A(°))*, = a, (A(°))^ = 0, (A(°))\ = 0, (A(°)r^.=/35i, (D2) 
with a > and /3 > 0. The square-root is expanded as 

VA = 5^ + e5(^) + 6^5(2) + 0(e3), (D3) 



where 



and 



(i?<°')', = (B'"')', =0, (Bl°')', = 0, (fiW)'^. = V^i', (D4) 



Here, i^^) ^ ^(2) _ 5(1) 
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